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The  Quiet  Sunrise  E  Region:  Enhancements  at  High  Latitudes 
in  Winter  Due  to  Increased  Nitric  Oxide 
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E  region  elcelron  eonecntriilUms  derived  from  Luropciin  Ineoherenl  Sealler  t  iicililv  obverviiiionv  ai 
70’N  were  found  to  he  much  larger  than  normal  for  high  volar  zenith  angicv  dining  .’A  of  A2  days  m 
winter.  Hleclron  concentrations  computed  with  the  Keneshea  code  are  compaiihle  with  the  ohsersa- 
tions  if  nitric  oxide  concenlrutions  are  about  lO'*  cm  *  at  HXH  10  kni  Lor  these  altitudes  near  ground 
sunrise,  it  is  shown  thtil  nunicricalty,  the  nitric  oxide  concentration  approximately  equals  the  sqiiaie 
of  the  electron  concentralion. 
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I.  iNIRtIUlK  ritlN 

Nitric  oxide  is  a  minor,  bm  important,  constituent  of  the 
atmosphere  up  to  approximately  150  km.  This  ubiquitous  gas 
helps  destroy  ozone  in  the  stratosphere  and  is  an  important 
IR  emitter  in  the  lower  thermosphere,  and  its  photoioniza¬ 
tion  creates  the  undisturbed  I)  region. 

In  the  region.  NO  has  iwci  clear  roles.  In  the  quiet 
nadday  E  region,  nitric  oxide  converts  O;'  ions  into  NO' 
ions  by  the  charge  transfer  process,  02  v  no  n() '  t 
Os.  This  reaction  increases  the  mean  recombination  coeffi¬ 
cient  of  the  E  region,  hut  only  slightly,  as  the  dissociative 
recombination  rate  coefficient  of  NO  '  is  about  twice  that  of 
Os'  and  as  the  electron  concentration  is  inversely  dependent 
upon  the  square  rout  of  this  coefficient,  rhesc  two  ions  are 
the  only  significant  ones  in  the  nonsporadic  E  layer.  By  day, 
they  have  comparable  abundances  [Kertcshco  el  til..  1970). 
At  night.  NO  '  is  by  far  the  more  dominant  ion.  arising 
principally  from  the  above  charge  transfer  process,  with  02 
produced  mainly  from  the  photoionization  of  Os  by  scattered 
solar  H  Ly  photons.  Some  NO  '  is  generated  by  the 
photoionization  of  NO  by  scattered  H  Ly  a  radiation.  (See 
Strahel  et  al.  [1980]  for  a  more  detailed  analysis  of  the 
various  faint  nighttime  ionization  sources,  including  star¬ 
light.) 

The  role  of  nitric  oxide  is  far  more  significant  in  the  E 
region  at  twilight  [SteiJer,  196.5;  .Swiilcr  and  Keneshea. 
I968|,  the  present  subject.  Near  ground  sunrise  (or  sunset), 
the  solar  H  Ly  «  flux  is  unattenuated  as  it  penetrates  into  the 
central  E  region,  I  lO-IOt)  km.  in  contrast  to  the  bulk  of  the 
daytime  radiation  (EUV  >-  X  rays)  which  ionizes O.  Ov.  and 
Nv.  Thus,  for  a  short  span  of  solar  zenith  angles  ^  near 
#  ground  sunrise  (or  sunset)  the  dominant  ionization  source  of 

the  E  regiixn  ts  NO  s  H  Ly  a  NO  ‘  +  e.  For  this  reason, 
the  ratio  of  NO  '  to  02  ions  is  greatest  at  twilight  {.Swkler 
*  and  Keneshea.  1%8;  Keneshea  et  al..  I97()|. 

2.  CaL(  lll.AtlON'.S 

rime-dependent  computations  of  electron  concentrations 
|<'l  were  accomplished  with  the  same  code  [Keneshea.  I%7| 
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used  for  the  first  time-dependent  model  of  the  h  region 
[Keneshea  el  a/..  I97()|.  Results  were  presented  (.'inider  el 
a/..  19761  but  were  ncxcr  submitted  previoush  for  publica¬ 
tion.  Neutral  concentrations  and  temperatures  were  taken 
I'rom  the  1966  U.S.  Standard  .Atmosphere  Supplements  for 
fall-spring  at  45  N  with  a  ItXKf'K  exospheric  temperature 
Calculations  were  performed  for  this  latitude  at  zero  solar 
declination  (equinox).  Our  NO  profile  is  based  mainly  on  (he 
observations  by  Meira  fI97lj  for  85-105  km.  with  extrapo¬ 
lation  to  higher  altitudes  by  diffusion,  and  henceforth  will  be 
denoted  as  the  baseline  profile.  )l  wilt  he  illustrated  later  in 
the  text  for  comparison  with  other  profiles,  fhe  NO  concen¬ 
trations  are  0.9.  1 .0.  and  0,9  in  units  of  lO'*  cm  '.at  IIX).  105. 
and  IK)  km.  respectively .  the  best  altitudes  for  comparing 
our  model  with  the  European  Incoherent  Scatter  Facility 
(EISCAT)  data.  Incidentally,  our  baseline  profile  closely 
matches  a  one-dimensional  model  for  NO  with  an  auroral 
electron  flux  of  0.5  ergs  cm  '  s  '  [Banh.  1992). 

•As  before  [Keneshea  et  al..  1970|.  scattered  intensities  for 
FI  Ly  <t  and  H  Ly  fi  at  night  arc  Kr  and  ().4'~r  of  their 
respective  noon  (.y  --  45°)  values.  Strahel  el  al.  ( 1974)  used 
more  detailed  scattering  profiles,  but  their  work  is  more 
pertinent  to  the  nighttime  E  layer.  These  scattered  intensi¬ 
ties  arc  of  little  significance  at  sunrise  but  only  establish  the 
nighttime  levels  of  the  electron  concentrations  prior  to 
sunrise.  .All  other  photon  intensities  also  itc  as  before, 
except  that  X  ray  intensities  were  reduced  by  a  factor  of  4. 
rhe  proper  level  of  X  ray  intensities  remains  somewhat 
controversial  (cf.  Siskind  et  al..  1990;  Clancy  et  al..  1992). 
This  problem  is  noi  of  concern  here,  as  X  rays  are  negligible 
in  the  lower  twilight  E  region  because  of  their  attenuation. 

Rate  coefficients  were  updated  from  those  listed  by  Ke¬ 
neshea  el  al..  [  1970).  and  in  particular,  the  rate  coefficient  for 
the  charge  transfer  process  02  t  NO  NO  ‘  +  0>  was 
lowered  slightly  from  6  to  4.6  x  10  "'em  '  s  However, 
since  the  ion  chemistry  of  the  lower  E  region  at  twilight  is 
reduced  basically  to  the  photoionization  of  NO  followed  by 
the  dissociative  recombination  of  NO  '  ions,  only  the  disso¬ 
ciative  recombination  coefficient  of  NO'  .  k.  is  significant. 
The  value  adopted  for  k  was  4  x  |0  '  cm  '  s  '  (.5(X)/T). 
Application  of  the  most  recent  value  for  the  rate  coefficient 
would  have  minimal  impact  on  the  calculations  as  it  differs 
little  from  that  used  here,  and  as  the  electron  concentration 
i.s  a  function  of  the  inverse  square  root  of  k. 

Production  rates  for  NO  '  at  twilight  from  solely  the 
photoionization  of  NO  by  H  Ly  a.  i/(NO.  H  Ly  n),  are 
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Fig.  1.  Productioa  rales  for  NO  '  at  twilight  from  the  photoioni/.Htion  of  NO  by  H  l.y  o  versus  (u  l  altitude  for  solar 
zenith  angles  in  2°  intervals  and  \h)  solar  zenith  angle  for  selected  altitudes.  iTen  limes  the  baseline  NO  concentrations 
used  in  the  calculations  are  given  in  Figure  .^.i 


displayed  in  two  formats  in  Figure  I .  Figure  1  h  shows  that  at 
1 10  km.  H  Ly  a  sunrise  begins  (sunset  ends)  near  y  =  94”. 
and  by  =  90°.  c/(NO.  H  Ly  a)  has  almost  attained  its 
maximum  value.  At  95  km.  H  Ly  a  sunrise  occurs  at  92.5°. 
and  <5f(NO.  H  Ly  a)  attains  nearly  its  maximum  value  by 
87.5°.  The  changes  in  <7! NO.  H  Ly  a)  are  quite  abrupt, 
almost  step  functions. 

The  ionization  coefficient  for  NO.  o-|/„.  is  6  x  10  ’  s  '. 
the  product  of  the  photoionization  cross  section,  ctj.  and  a 
nominal  photon  flux  of  3  x  10 "  cm  '  s^'  for  /o.  the 
intensity  of  H  Ly  a.  The  maximum  ionization  rate  of  NO  by 
H  Ly  a  is  nr[/|,[NOJ.  At  105  km.  for  example,  the  product  is 
60  cm  s  '.  This  NO*  production  rate  is  a  factor  of  5 
higher  than  that  calculated  by  Swider  and  Keneshea  ( 1%81, 
since  they  used  values  of  3.2  x  lo’  cm  for  [NO]  and  2  x 
10"  cm'  ^  s“'  for  /,).  or  a  value  of  0.21  lower  for  /,|[N01. 
This  lower  value  for  /,j  is  more  appropriate  for  solar  mini¬ 
mum. 

3.  CoMPAKi.soNS  With  Expkrimkntai.  Data 

The  EISCAT  radar  has  provided  evidence  of  an  enhanced 
E  layer  during  quiet  conditions  for  x  <  93°  [Kirkwood.  1991  j. 
The  layer  was  present  on  809F  of  32  winter  days.  Electron 
concentrations  \e\  from  F3ISCAT observations  are  compared 
in  Figure  2  with  |e)  profiles  computed  with  our  baseline  NO 
profile  (Figure  2o)  and  with  NO  concentrations  lOX  (Figure 
2b)  and  l(K)X  (Figure  2< )  greater.  From  comparison  of 
Figure  la  with  Figure  I  a.  the  other  ionization  sources  (EUV 
+  X  rays)  become  important  above  1 10  km  as  x  decreases 
below  90”.  For  Figure  2c,  ionization  of  NO  is  the  dominant 
ionization  mechanism  at  all  altitudes.  Note  that  near  100-1 10 
km.  computed  [<'|  increases  by  a  factor  of  3  for  an  [NO) 
increase  of  a  factor  of  10  at  altitudes  where  <y(NO,  H  Ly  a) 
is  the  dominant  ionization  source,  since,  as  shown  later. 


[NO)  is  roughly  proportional  to  jclv  The  NO  concentra¬ 
tions  appropriate  to  the  EISCAT  measurements  arc  esti¬ 
mated  from  Figures  Ih  and  2<  to  reach  concentrations 
slightly  above  It)'*  cm  '  at  l(K)-l  10  km.  i.e.,  about  a  factor 
1.5-30  limes  greater  than  our  baseline  levels  for  nitric  oxide. 
For  8.5-95  km,  the  increase  in  |NO|  is  a  factor  of  10  or  less 
greater  than  our  baseliiie  nitric  oxide  concentrations. 

The  functional  relationships  for  NO  arc  better  understood 
if  we  consider  a  simple,  steady  state  solution,  permissible 
because  of  the  slow  rate  of  change  of  the  solar  zenith  angle. 
l,5°/min,  appropriate  to  the  EISCAT  data:  70°N;  mid- 
November;  X  ~  87”.  The  steady  state  solution  for  NO  is 


where  t  =  [OyJ  H  Ch(y.  X).  By  happenstance,  this 
equation  simplifies  to 

(NO)  =  [(')”  exp  (t)  (2) 

since  the  ionization  coefficient  for  NO,  fr,/,,.  is  6  X  10  s  '. 
as  noted  above,  and  since  A  -  6  x  10  ’  cm  '  s  '.  T 
averaging  about  2I0°K.  The  absorption  cross  section  of  O;  at 
121,6  nm  (H  Ly  «),  rr^,  is  10  cm  N  The  average  scale 
height  H  at  8.5-110  km  is  5.8  km.  The  Chapman  function. 
Ch(87°.  A"),  is  about  15.6  for  X  =  (,■;-*  R)iH.  where  .-(mean) 
=  95  km  and  /?( Earth)  =  6370  km. 

Using  the  EISCAT  results  for  [(>).  equation  (2)  was  solved 
for  [NOl  at  85.  90,  95,  100.  105.  and  1 10  km.  The  solutions, 
lower  bounds  compared  to  the  time-dependent  solutions,  arc 
plotted  as  a  smooth  line  labeled  87°  in  Figure  3  in  conjunction 
with  our  baseline  [NO]  distribution.  Lower  |N()|  is  calcu¬ 
lated  if  T  and/or  /„  arc  greater,  according  to  equation  (I). 
Different  models  of  the  atmosphere  mainly  would  shift  the 
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Fig.  2.  Electron  concentrations  (<>]  versus  altitude  obtained  from  EISCAT  observations  (bold  curve)  at  x  =  b7' 
[Kirkwood  t';<9tl,  compared  with  computations  of  (e]  performed  with  various  NO  concentrations;  (a)  baseline  |N01 
and  with  [NO]  increased  by  factors  of  [h)  10  and  (<  )  100  limes  the  baseline  values.  (Ten  times  the  baseline 
concentrations  are  plotted  in  Figure  3.) 


entire  computation  up  or  down  a  km  or  two.  It  should  be 
particularly  noted  that  for  the  al'itudes  100-110  km.  the 
solution  is  trivial,  as  the  attenuation  of  H  Ly  o  is  negligible. 
Equation  (2)  is  reduced  to  merely 

(NOJ  =  [e]’  (3) 

Kirkwood  [\99\  \  calculated  an  electron  concentration  of  4 
X  10-’  cm'^  for  100  km  at  ;^  =  87°  for  (NO)  =  9.5  x  10^ 
cm a  value  taken  from  Win  et  al.  (1976).  whose  (NO]  for 
100-110  km  are  within  20%  of  our  baseline  concentrations. 
She  calculated  (e)  =  9  x  lO^  cm'^  by  increasing  Witt  et  al's 
NO  concentration  by  a  factor  of  5.  The  increase  from  4  to  9 


Fig.  3.  Nitric  oxide  concentrations  (bold  curve)  estimated  from 
EISCAT  electron  concentrations  at  =  87°  using  equation  (2).  The 
thin  curve  is  the  baseline  [NO]  profile  used  in  our  computations 
increased  by  a  factor  of  10, 


is  in  line  with  the  simple  theory  above  (equation  (3)1. 
However,  we  calculate  (e)  =  9  x  10 '  cn  '  for  our  baseline 
NO  profile  (Figure  2a).  suggesting  that  her  computations 
involve  excessive  attenuation  of  H  Ly  o.  We  estimate  that 
Kirkwood's  [1991)  computations  would  require  NO  values 
near  lO'®  cm'^  to  match  the  EISCAT  data  at  !(X)-1I0  km. 

4.  Discussion  and  Conclusions 

We  proposed  over  a  decade  ago  \Swider  et  at..  1976.  1977] 
that  analysis  of  the  high-latitude  quiet  E  region  at  sunrise 
might  yield  enhanced  concentrations  of  nitric  oxide.  The 
EISCAT  results  reported  by  Kirkwood  [1991)  appear  to 
support  this  prediction.  Both  our  lime-dependent  and  rough 
steady  state  calculations  indicate  that  NO  concentrations 
reach  1-3  x  lO"^  cm"*  at  100-110  km.  This  enhancement 
provides  for  an  augmented  E  region  at  sunrise  since  H  Ly  a 
photons  essentially  penetrate  these  altitudes  without  any 
attenuation  for  solar  zenith  angles  near  90°  [Swider  and 
Keneshea,  1968).  We  have  shown  that  near  ground  sunrise, 
;y  =  90°  the  nitric  oxide  concentration  approximately  equals 
(numerically)  the  square  of  the  electron  concentration  at 
lOO-l  10  km.  This  relationship  is  true  for  all  latitudes  under 
quiet  conditions,  providing  the  nitric  oxide  concentration  is 
not  too  low,  i.e..  that  H  Ly  a  +  NO  -♦  NO*  +  e  is  the 
dominant  ionization  mechanism. 

In  the  polar  regions,  all  other  factors  being  equal.  NO 
concentrations  are  greatest  in  winter,  when  photodissocia¬ 
tion  of  NO  is  inoperative.  Nitric  oxide  concentrations  de¬ 
rived  from  satellite  -yband  observations  |/?w.v(  6  and  Barth, 
1975)  show  (NO)  enhanced  in  the  polar  regions  during 
winter.  The  main  surprise  in  the  EISCAT  observations  may 
be  that  the  (NO)  enhancements  were  present  72%  of  the 
time,  23  of  32  cases  [Kirkwood.  1991).  Nitric  oxide  concen¬ 
trations  derived  from  .satellite  measurements  of  the  NO  y 
band  [e.g.,  Feaen  et  al..  1990;  Barth.  1992)  appear  to  be 
tower,  on  average,  than  those  estimated  here  and  those 
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derived  from  in  situ  aurora!  ion  composition  measurements 
[Sw’ider  anii  Niinisi,  1977],  Although,  in  the  latter  study,  in 
at  least  three  of  the  eight  auroral  events  considered.  fNO| 
equaled  10*^  cm  '  at  100  km.  In  part,  the  lower  values  for  the 
satellite-based  data  may  be  attributable  to  their  being  aver¬ 
aged  over  a  much  wider  area  than  either  the  EISCAT-  or  the 
rocket-based  data. 

We  have  not  presented  any  sunset  profiles  because  their 
interpretation  is  more  comple.x.  Perturbations  by  sporadic  E 
and  aurora  are  more  difficult  to  discriminate  from  the  normal 
daytime  E  region  behavior  in  comparison  to  the  low'  back¬ 
ground  conditions  prior  to  sunrise.  However,  the  late  after¬ 
noon  data  from  EISCAT  are  quite  similar  lKirh<.<H>d.  1991) 
to  the  early  morning  data  depicted  here. 

Other  mechanisms  for  an  enhanced  E  region  at  sunrise  are 
thought  less  likely.  There  is  no  evidence,  for  e.xample.  of 
enhanced  layers  of  meteoric  ions  at  high  latitudes,  although 
the  data  ba.se  is  meager  [Slider,  1984], 

AcknowledRment.  The  Editor  thanks  C.  G,  Fesen  nd  another 
referee  for  their  assistance  in  evaluating  this  paper. 
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